Ab initio study of spin-dependent transport in carbon nanotubes with iron and 

vanadium adatoms 
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We present an ab initio study of spin dependent transport in armchair carbon nanotubes with 
transition metal adsorbates, iron or vanadium. We neglect the effect of tube curvature and model 
the nanotube by graphene with periodic boundary conditions. A density functional theory based 
nonequilibrium Green's function method is used to compute the electronic structure and zero-bias 
conductance. The presence of the adsorbate causes a strong scattering of electrons of one spin 
type only. The scattering is shown to be due to coupling of the two armchair band states to 
the metal 3d orbitals with matching symmetry causing Fano resonances appearing as dips in the 
transmission function. The spin type (majority/minority) being scattered depends on the adsorbate 
and is explained in terms of d-state filling. The results are qualitatively reproduced using a simple 
tight-binding model, which is then used to investigate the dependence of the transmission on the 
nanotube width. We find a decrease in the width of the transmission dip as the tube-size increases. 

PACS numbers: 73.63.Fg, 73.63.-b 
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I. INTRODUCTION 

Various carbon structures have attracted considerable 
attention in the last decade or so, due to their po- 
tential use within future nanodevicesi. Graphene, the 
single-atom thick two-dimensional sheet of graphite is 
the basis for many carbon materials. Rolling up such 
a sheet creates the one-dimensional single-walled carbon 
nanotube (SWNT) and wrapping up a sheet makes OD 
fuUcrenes. These materials have a wide range of remark- 
able electronic properties. Recently graphene ribbons - 
cut graphene sheets - were used to make the first tran- 
sistor carved in graphene^ which is stable at room tem- 
perature. 

Not only pushing the limits within electronics, these 
carbon materials may lay the ground for the emerging 
field of spintronics^i^, where the functionality of the 
device is based on the spin degree of freedom. While 
graphene and SWNT's are not inherently magnetic, in- 
troducing defects^iii^, impurities or boundaries^iii^ may 
alter this fact. In recent years several studies have 
been carried out on SWNT/graphene - transition metal 
adatom systems. First principles calculations have 
been employed to determine equilibrium geometries and 
magnetic properties: Fagan et al. reported electronic 
structure studies on single iro n^^d^d^ and manganese 
adatomsiiii^ on zigzag SWNT's showing a total magnetic 
moment of the systems and in most cases a magnetisation 
of the tube itself. A similar earlier study with additional 
transition metal elements on graphene was published by 
Duffy et ali^ reaching similar conclusions. Manganese 
dimers, trimers and wires on zigzag SWNT's^^ exhibit 
magnetic moments close to free manganese. Yang et 
alJ^ have proposed that iron and cobalt coated and filled 



SWNT's can work as spintronic devices demonstrating a 
spin polarisation close to 90 % at the Fermi level as well 
as considerable magnetic moments. Interestingly, both 
semiconducting and metallic SNWT's showed high spin 
polarisation. Kang et aZ^— obtained similar conclusions 
for various iron nanowire configurations inside armchair 
nanotubesi^. On the other hand, work by Kishi et a/.— 
shows that wires of iron and cobalt lose their magnetic 
moments when adsorbed on armchair SNWT's. Filling of 
SWNT's with transition metals has been realised experi- 
mentally several years agoi^, which adds to the potential 
of these systems. A few studies have been published on 
electron transport in this context. Iron-SWNT's junc- 
tions with Ceo molecule&lS and pristine SWNT's^" have 
been proposed as magnetic tunnel junctions. First prin- 
ciples transport calculations yield tunnel magnetoresis- 
tances of 11 % and 40 %, respectively, for these systems. 
A spin-polarised current has also been reported in theo- 
retical works on nanoribbons with substitutional boron 
atoms^^, where spin-dependent scattering is found. 

Motivated by the reported spin-dependent scattering, 
we seek here to explain the physical mechanisms found in 
the numerical calculations. We investigate the influence 
on the transport properties of a single iron or vanadium 
atom adsorbed on a SWNT. By performing ab initio spin- 
polarised transport calculations we demonstrate a spin- 
dependent scattering. The spin type being scattered is 
shown to depend on the type of adsorbate. By analysing 
the PDOS of the 3d orbitals of the adsorbates we find the 
scattering to be caused by coupling of the band states in 
the SWNT to these orbitals, resulting in Fano resonance 
phenomena. Using a simple tight-binding model we qual- 
itatively reproduce these results, and also investigate the 
dependence of the transmission on tube size. 
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FIG. 1: (Colour online) The unit cell of the system. The 
sheet is cut in an armchair structure along the y-direction. 
The adatom, iron or vanadium is placed in the center. The 
shaded areas marked L and R are the electrodes, and C is the 
central region. Periodic boundary conditions are imposed in 
the transverse direction (x). 



The paper is organised as follows. Sec. |TT] introduces 
our model system and the technical details of the calcula- 
tions. The ab initio results are presented in Sec. Illll and 
the Fano resonance is briefly introduced along with trans- 
mission eigenchannel analysis. In Sec. IIVI we present a 
simple tight-binding model which is compared to the ah 
initio results. Results for larger tubes are given based on 
the simple model. 



II. SYSTEM 

The system used for transport calculations is shown 
in Fig. [T] The graphene sheet is cut in an armchair 
structure along the y-direction which results in a metal- 
lic system. We employ periodic boundary conditions in 
the transverse x-direction. We have used the ab ini- 
tio pseudopotential density functional theory (DFT) as 
implemented in the SIESTA code^^ to obtain the elec- 
tronic structure and relaxed atomic positions from spin- 
polarised DFT. We employ the GGA PBE pseudopoten- 
tial for exchange-correlation^^. Our spin transport cal- 
culation is based on the nonequilibrium Greens function 
method as implemented in the TranSIESTA^i code, ex- 
tended to spin-polarised systems^^. However, we only 
consider here the zero bias limit and focus on electron 
transmission close to the Fermi energy. According to the 
TranSIESTA method^! the system is divided into left and 
right electrodes, marked L and R on Fig. [1] and a central 
region marked C. A single iron or vanadium adatom is 
placed in the middle of the central region. The electrodes 
both contain 32 atoms while the central region consists 
of 64 sheet atoms. The adatom systems are relaxed using 
the CG method with a force tolerance of 0.01 eV/A. The 
carbon atoms were kept fixed during geometry optimisa- 
tion. The relaxed position of the adatoms is in the center 
of a hexagon at a distance of 1.73 A and 1.89 A from the 
sheet plane for iron and vanadium, respectively. The pe- 
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FIG. 2: (Colour online) The transmission of the vanadium 
adatom system using a single-zeta basis set for the whole sys- 
tem versus using double-zeta polarised for the adatom. 



riod of the cell in the transverse x direction is 8.52 A, 
which is then the smallest distance between adatoms 
in neighbouring cells, and should be enough to prevent 
significant adatom - adatom interactions. We employ 
Gamma point sampling (fc = 0) in the periodic x direc- 
tion. The system then corresponds to an (2,2) armchair 
nanotube described by the approximate Graphene Sheet 
ModeP^, that is, neglecting curvature effects. We will re- 
turn to the dependence of the electron transmission close 
to the Fermi energy on the GSM-armchair tube width 
(n, n) in Sec. IIVI The mesh cutoff value defining the 
real space grid was set to 175 Ry. A single-zeta (SZ) ba- 
sis set was used. The difference between a SZ basis for 
all atoms and increasing the vanadium basis to double- 
zeta polarised (DZP) can be seen on Fig. [21 The only 
qualitative difference between using DZP and SZ is the 
additional dip (a single point in the graph) at i? ~ 0.25 
eV in the SZ case. The reason for this not occuring for 
DZP is merely a matter of resolution. The dips will be 
addressed in the following two sections. In the remains 
of the paper we have used a single-zeta basis. 



III. RESULTS 

The spin resolved transmissions as a function of energy 
relative to the Fermi level, Ep, of the graphene sheet 
with iron and vanadium adatoms are shown in Fig. [3l 
The transmission of a pure sheet is 2 for each spin type, 
since there are two bands in the energy window each con- 
tributing with a fully transmitting channel for each spin. 
It is seen that spin dependent scattering occurs due to 
the presence of the adatoms. In the case of iron the mi- 
nority spin type is significantly suppressed around Ep^ 
whereas the majority spin electrons transmit completely. 




FIG. 3: Top: Transmission as a function of energy for iron (left) or vanadium (right). Within a 0.5 eV range of Ep the minority 
spin channels are suppressed in the case of iron whereas the majority channels are suppressed for vanadium. Middle: The 
transmission of the two minority (majority) spin channels for the iron (vanadium) system. The sum of the two channels yield 
to total transmission in each case. Bottom: The projected density of states (PDOS) of the 3d orbitals of the iron(vanadium) 
adatom for minority (majority) spin. 



For vanadium, likewise, scattering occurs for only one 
spin channel, but in this case it is the majority spin elec- 
trons which are scattered. The transmission of the two 
bands for minority (iron) and majority (vanadium) spin 
is shown on Fig. [3l middle graph. We see that each band 
closes completely at certain energies. 

In the case of vanadium {ScP-is^) we expect from 
Hund's rules a total spin of 3 and the majority d-states 
will be located around Ep^ whereas the minority states 
are all empty and well above Ep. The MuUiken analysis 
indicates a half filling of the majority spin 3dyz and Sd^x 
orbitals and a full filling of Sd^2_y2, "idxy and 'id^z2_j.2. 
For minority electrons all id orbitals are empty. This is in 
full correspondence with the 3(i orbital PDOS plot shown 
in Fig. [3] bottom graph. In the case of iron (3d®4s^) 
we expect from Hund's rules that the majority states 
are all filled and well below Ep, and now the partially 
filled minority states are located around Ep yielding a 



total spin of 4. This is again supported by the MuUiken 
analysis data as well as the projected density of states 
(PDOS). Comparing PDOS and the transmission there 
is a strong correlation between adatom orbital energies 
and conductance dips. This is very clear for the two sep- 
arated channel closings for vanadium a,t E k, Ep and 
E — Ep ss —0.5 eV. The interference between a waves 
involving the quasi-bound d-state on the adatom and di- 
rectly transmitted band states yield a Fano (anti-) reso- 
nance. The line shape of the resonance is given by the 
Fano functional 

m = (1) 

where e = {E—Eb) /V . Here F is the resonance width, Epi 
the resonance energy and q the asymmetry parameter. 
Despite the overlapping of dips in our calculations their 
shape appears symmetric in the vicinity of zero transmis- 
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FIG. 4: (Colour online) The real part of the (left-to-right) eigenchannel scattering states for the vanadium system, stemming 
from majority spin electrons from both bands, at the energy of the first dip (-0.5 eV) (a),(b) and second dip (-0.1 eV) (c),(d) 
in the transmission spectrum. Only the solutions in the scattering region are shown. The size of the shapes indicates a cut-off 
value for the wavefunction. White indicates a positive sign a blue a negative sign of the wavefunction. The involved orbitals 
of the vanadium atom (marked with dotted circle) are seen to be (a) : 3d^2_y2 (and a minor presence of 3(1-^^2 _r2), (b) : 3dxy, 
(c): 3dyz and (d) : 3dzx- 



sion. Since the structure has inversion symmetry around 
the adatom in the transport direction the parameter q 
must be real and thus we expect to have zero transmission 
dips as seen in Fig. [3] at the anti- resonances e = —q^i^. 
A very sharp resonance can be seen for vanadium around 
— leV which should also go to zero. However, it is re- 
lated to the most localized d-orbital, c?3^2_^2, and a small 
asymmetry in the numerical calculation and the resolu- 
tion cause it not to go exactly to zero. Signatures of Fano 
resonances have previously been observed experimentally 
in the conductance of multi-wall carbon nanotubes at low 
temperature^. 

The anti-resonances are illustrated further by per- 
forming an eigenchannel analysis where scattering states 
with well-defined transmissions (" eigenchannels" ) are 
constructed^. We take vanadium as our example. 
The eigenchannel transmissions, T'„(e), provide the 
transmission of channel n at a given energy e. Plotted 
on Fig. m (a)(b) and (c)(d) are the real part of the 
majority spin scattering state solutions of both bands 
for energies corresponding to the first dip (— 0.5eV) and 
the second dip (— O.leV), respectively. We see that the 
wavefunction is nonvanishing only on the left side of the 
vanadium atom as expected since we have full reflection 
in all cases. The anti-symmetric and symmetric solutions 
are seen on Fig. [11(a) (d) and (b)(c), respectively. Note 
that in the transverse direction the neighbouring C 



atoms have the same (opposite) sign in the symmetric 
(anti-symmetric) case. These solutions are indeed 
matched by the vanadium atom orbitals, which is seen 
on the figures by a match of signs (colours). From the 
shape and sign in the plots the orbitals involved can be 
identified for each band at both energies. 



IV. A TIGHT-BINDING MODEL 

We will now rationalize the results for the channels 
and transmissions in terms of the simplest possible tight- 
binding model. We consider only the coupling of the 
adatom rf-orbitals with the 6 nearest 7r-orbitals. The 
starting point will be the two band states of the tt- 
electrons in the armchair direction at the Fermi level, see 
Fig. [51 These are characterized by rotational symmetry 
around the tube axis and come in an odd/even version 
around the symmetry plane normal to x along the tube 
(y) and cutting through the adatom. They couple to 
different d-orbitals on the adatom with the same sym- 
metry. Thus the symmetric/anti-symmetric band only 
couple to the d-orbitals even/odd in x. To make a min- 
imal model (with a minimum of parameters) we assume 
that the d-orbitals have the same on-site energy Ed, and 
a coupling to the carbon 7r-orbitals described in units of 
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FIG. 5: (Colour online) Top panel: The two vr-band (pz) 
states (not normalized) in a hexagon where the arrow de- 
notes the direction along the armchair tube {y direction) with 
= e'"'^^^^. Left/right panel correspond to symmetric/anti- 
symmetric solution on the A,B-dimers {x direction). Lower 
panels: Simple tight-binding model. The sign and relative size 
of the hopping matrix elements between the ring 7r-orbitals 
and the metal adatom d-orbitals when the metal-atom is sit- 
uated I.7A above the middle of the ring. 



Vpd = Vpdo- taking Vpdrr « O.bVpdcr^- The size and sign of 
the coupling of the different d-orbitals to the Pz orbitals 
are illustrated in Fig. [5l left/right panels corresponding 
to coupHng to the two types of bands. For both bands 
it is seen that the d^z^-r^ orbitals do not couple at Ep 
since the wavefunction values around the hexagon sums 
up to zero. Therefore we do not expect significant con- 
tributions from this orbital, which will be very localised. 

We only consider the majority spin d-statcs for V, since 
they, as mentioned above, will be located around Ep with 
Ed ~ 0. Likewise, for iron we only consider the minority 
spin states. The model calculation with {Ed, Vpd) pa- 
rameters chosen to fit the data from the full ab initio 
calculation arc shown in Fig. [SI It is seen that the trans- 
mission is more suppressed in the full calculation where 
the individual d-orbitals are allowed to have different en- 
ergies, but the fact that we see 4 dips corresponding to 
the four coupling d-orbitals is clear. 

With the simple model we can now explore what hap- 
pens when the diameter of (n, n)-tube is increased. This 
is shown in Fig. [7] where it is seen that the transmis- 
sion still goes to zero at the anti-resonance points but 
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FIG. 6: (Colour online) The transmission through the 
symmetric/anti-symmetric bands for Fe minority spin (left 
panels) and V majority spin (right panels) calculated with 
the simple tight-binding model. Lower panels display the pro- 
jected density of states onto the d-orbitals. 



in an interval getting narrower with increasing n. The 
latter is in agreement with the findings of Todorov and 
White^. They conclude that the larger the tube diame- 
ter the more the scattering due to an impurity is reduced. 
This happens because the overlap with the wavefunc- 
tion (normalised around the tube circumference) with 
the scattering potential goes down. The fact that the 
conductance still drops to zero in our calculations is due 
to the Fano resonance phenomena. 



V. CONCLUSION 

We have described spin-polarised zero bias transport 
calculations for armchair carbon nanotubes with ad- 
sorbed single iron or vanadium atoms. We find a signif- 
icant difference between transmissions of majority and 
minority spin. The presence of the metal adatoms causes 
spin-dependent closing of the conduction channels at cer- 
tain energies. The mechanism is due to Fano resonances 
related to the particular d-states close to the Fermi en- 
ergy. Only d-orbitals with a symmetry matching the 
symmetry of the Bloch band solutions take part in the 
scattering. The scattered spin type (minority or major- 
ity spin) can be explained by the filling of the 3d orbitals 
via Hund's rule. Results of a simple tight-binding model 
show that increasing the width of the tube, and thus re- 
ducing the concentration of adatoms, results in narrower 
dips in the conduction, but still with complete closing. 
The latter is explained by the reflection symmetry of the 
system in the transport direction. 
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FIG. 7; (Colour online) The dependence on tube-size, (n, n) 
calculated with the simple tight-binding model at the first dip 
in the Fe case. 
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